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We discuss the ground state properties of a spin 1/2 magnetic ion with threefold t2g orbital 
degeneracy on a highly frustrated pyrochlore lattice, like Ti'^^ ion in B-spinel MgTi2 04. We for- 
mulate an effective spin-orbital Hamiltonian and study its low energy sector by constructing several 
exact-eigenstates in the limit of vanishing Hund's coupling. We find that orbital degrees of freedom 
modulate the spin-exchange energies, release the infinite spin-degeneracy of pyrochlore structure, 
and drive the system to a non-magnetic spin-singlet manifold. The latter is a collection of spin-singlet 
dimers and is, however, highly degenerate with respect of dimer orientations. This "orientational" 
degeneracy is then lifted by a magneto-elastic interaction that optimizes the previous energy gain 
by distorting the bonds in suitable directions and leading to a tetragonal phase. In this way a va- 
lence bond crystal state is formed, through the condensation of dimers along helical chains running 
around the tetragonal c-axis, as actually observed in MgTi2 04. The orbitally ordered pattern in the 
dimerized phase is predicted to be of ferro-type along the helices and of antiferro-type between them. 
Finally, through analytical considerations as well as numerical ab-initio simulations, we predict a 
possible experimental tool for the observation of such an orbital ordering, through resonant x-ray 
scattering. 



PACS numbers: 75.10.Jm, 75.30.Et 
I. INTRODUCTION. 



A frustrated antiferromagnet is characterized by the 
topology of underlying lattice and/or by the presence of 
competing interactions that preclude every pairwise mag- 
netic interaction to be satisfied at the same time. Such 
physical systems have recently attracted a wide interest 
due to the concept of "macroscopic" degeneracy, namely 
the existence of a huge number of states with the same 
energy. This degeneracy in the ground state manifold 
can be usually removed through a large variety of ef- 
fects, like order-out-of-disorder mechanisms^, by thermal 
or quantum fluctuations, or spin-Peierls like symmetry- 
lowering transitions."^'^ In some highly frustrated mod- 
els, like nearest-neighbor Heisenberg antiferromagnet on 
pyrochlore lattice, the former mechanisms are inactive^ 
and such a spin system would remain liquid down to the 
lowest temperatures^^ unless magnetoelastic couplings in- 
duce a symmetry-breaking transition. Yet, in real com- 
pounds geometrical frustration can also be partially or 
fully released when magnetic ions forming a frustrated 
lattice possess an orbital degeneracy. Typically this hap- 
pens in transition metal ions with orbitally degenerate 
partly-filled d-levels. The physical behavior of such sys- 
tems is expected to be drastically different from that of 
pure spin models, as the occurrence of an orbital ordering 
(00) can modulate the spin exchange and lift the geo- 
metrical degeneracy of the underlying lattice. Indeed, a 
variety of novel phases driven by orbital degrees of free- 
dom can be stabilized in this way. Among known exam- 



ples are vanadium cP compounds with frustrated lattices, 
where the orbital order is shown to induce a spin-singlet 
ground state without any long-range magnetic order for 
triangular lattice ^2. or a spin ordered one, for pyrochlore 
latticeiS 

In this paper we study a system with threefold- 
orbitally-degenerate S = 1/2 magnetic ions in a corner- 
sharing tetrahedral (pyrochlore) lattice. Our work is 
motivated by the very recent synthesisSi and interesting 
experimental data on iJ-spinel MgTi204f2ii2i a d-'^-type 
transition metal oxide. Here magnetically active Ti"^"*" 
ions form a pyrochlore lattice and are characterized by a 
single electron in a c?-shell. The crystal field of the oxy- 
gen octahedron surrounding each Ti ion splits this d-level 
into a high energy doublet eg and a low energy triplet t2g 
in which c?-electron residesi^i The ground state of a spin 
one-half Ti'^^ ion is thus threefold-orbitally-degenerate. 
The recent experiments on MgTi204 have shown that 
this compound undergoes a metal-to-insulator transition 
on cooling below 260 K, with an associated cubic-to- 
tetragonal lowering of the symmetry.^ At the transition 
the magnetic susceptibility continuously decreases and 
saturates, in the insulating phase, to a value which is 
anomalously small for spin 1/2 local moments: for this 
reason the insulating phase has been interpreted as a 
spin-singlet phase. Subsequent synchrotron and neutron 
powder diffraction experiments have revealed that the 
low-temperature crystal structure is made of alternating 
short and long Ti-Ti bonds forming a helix about the 
tetragonal c-axisiiS, These findings have suggested a re- 



2 



moval of the pyrochlore degeneracy by a one-dimensional 
(ID) helical dimerization of the spin pattern, with spin- 
singlets (dimers) located at short bonds. This phase can 
be regarded as a valence bond crystal (VBC) since the 
long-range order of spin-singlets extends throughout the 
whole pyrochlore lattice. 

The aim of the present work is to discuss the micro- 
scopic mechanism behind the realization of this unusual 
and intriguing VBC structure on the pyrochlore lattice. 
We argue that the key role in this mechanism is played by 
orbital degeneracy and, remarkably, one does not have to 
invoke any additional exotic interaction to stabilize such 
a novel phase, as necessary for purely spin models. In- 
deed, due to the orbital degeneracy, various type of spin- 
singlet phases, such as resonating valence bond (RVB) 
and VBC states can be formed even for unfrustrated cu- 
bic latticeii^ 

The orbital degree of freedom does modulate spin ex- 
change energies, thus removing the infinite spin degen- 
eracy, characteristic of pyrochlore structures, and drives 
the system to a non-magnetic spin-singlet state. This lat- 
ter is a collections of spin-singlet dimers with a residual 
macroscopic degeneracy of orientational character. The 
residual degeneracy is then lifted by a magnetoelastic in- 
teraction, that optimizes superexchange energy gain by 
distorting each tetrahedron in such a way as to lead to the 
experimentally observed helical pattern in MgTi204. We 
also find that the helical dimerized state is accompanied 
by a peculiar orbital pattern in which orbital order is of 
ferro-type along the helices and of antiferro-type between 
them. 

Moreover, we can show how to identify such an orbital 
ordering experimentally, by means of resonant x-ray scat- 
tering (RXS). In fact, differently from what happens in 
the case of manganitesfi^ti^ii^ii^ where the ratio between 
the 00-induced and Jahn- Teller (JT) induced effects is 
about 1/10 in amplitude, for MgTi204 this ratio is about 
1/3, due to the less distorted oxygen octahedra, as well 
as to the presence of t2g electrons instead of Cg, that 
couple less to the oxygen environment. Such a reduced 
ratio allows a subtle interference effect between 00 and 
JT terms, that gives rise to an increase of the signal by a 
factor of about 1.7. We have performed a detailed analyt- 
ical analysis as well as an ab-initio numerical simulation 
to suggest some experiments in this direction. 

In more detail, the paper is organized as follows: in 
Section II we derive an effective Kugel-Khomskiipii model 
Hamiltonian, and deduce its possible low-energy states 
on the pyrochlore lattice. We then introduce the spin- 
singlet ground state manifold and discuss its degeneracy 
due to the dimer coverings of the pyrochlore lattice. In 
Section III we consider the magnetoelastic interactions 
of the pyrochlore lattice and analyze the effects of the 
coupling of bond distortions with the orbitally driven ex- 
change modulation, first qualitatively, and then quantita- 
tively, underlining the infiuence of the spin-singlet phase 
on the strength of such a distortion. In Section IV the ef- 
fect of an applied magnetic field is briefly discussed, and. 



finally, in Section V we analyze the crystal structure and 
the orbital symmetry of the tetragonal phase, deriving 
the analytic expression for the measurable quantities in 
a RXS experiment. We also perform a series of numeri- 
cal ab-initio simulations by means of the finite difference 
method, implemented in the FDMNES packagCfi^ in or- 
der to propose a possible experiment to detect the orbital 
pattern. To faciliate the reader, two Appendices, A and 
B, are given for technical details. 

Part of the results presented here were already an- 
nounced in a previous short communication^ 



II. MODEL AND FORMALISM. 
A. Effective spin-orbital Hamiltonian. 

Here we discuss the superexchange spin-orbital Hamil- 
tonian for threefold orbitally-degenerate d^-ions on a py- 
rochlore lattice. We assume that the insulating phase 
of MgTi204 is of Mott-Hubbard type and can, thus, be 
described by the Kugel-Khomskii modeliii We consider 
the system in its cubic structure and look for possible 
instabilities towards symmetry reductions. 

The relevant electronic degrees of freedom are de- 
scribed by spin S = 1/2 and pseudospin r = 1 op- 
erators. This latter labels the orbital occupancies of 
t2g orbitals, (|a/?) = l^;^), \xz), \yz)), with the cor- 
respondence: = — 1 \yz), = — > \xy), and 
= 1 — !■ \xz). Our parameters are the nearest-neighbor 
(NN) electron hopping matrix t, defined in the Appendix 
A Eq. (|A1I) . the Coulomb on-site repulsions Ui (within 
the same orbital) and U2 (among different orbitals), and 
the Hund's exchange, Jh- For t2g wavefunctions the re- 
lation Ui = U2 + 2Jh holds due to rotational symmetry 
in real space. 

Considering the pyrochlore structure the effective 
Hamiltonian can be simplified by retaining only the lead- 
ing hopping parameter, i, due to the NN dda overlap and 
neglecting the smaller ddir and ddS contributions. This is 
justified by the fact that the transfer integrals due to the 
TT and d bonding are, respectively, around 1/10 and 1/3 
of that of a bonding. S2, The major advantage of this sim- 
plification, is that dda overlap in a/3 plane connects only 
the corresponding orbitals of the same aP type. Thus, 
the total number of electrons in each orbital state is a 
conserved quantity and, therefore, the orbital part of the 
effective Hamiltonian (|A2p becomes Ising-like (or better, 
a Potts Zs-like), making possible to get analytical results 
about the ground-state. The spin-orbital Hamiltonian is 
presented in Appendix A in the most general form [see 
Ea. (|A2|l ]. when also ddn and dd6 hopping elements are 
considered. Here we introduce the simplified version with 
only dda hopping terms: 

i?cff = -JiJ2[^^- ^1 + 3/4] 



+J2 [S^ ■ Sj - 1/4] + J3 ^ [s^ ■ S, - 1/A]d,, 

where the sum is restricted to the NN sites on the py- 
rochlore lattice. The orbital contributions along the bond 
ij in a/?-plane is given by 

Oij ~ Pi,al3ii — Pj,a(}) + ^j-, 0/3(1 — Pi^afi) 

Oij — Pi,al3Pj,al3, (2) 

where Pi^ap — T"-"- stands for the projector on orbital 
state I a) — \aP) and is defined in Appendix A. 

The first and second terms in Hcs ^ describe the 
ferromagnetic (FM) Ji — /{U2 — Jh) and the anti- 
ferromagnetic (AFM) J2 — t^/{U2 + Jh) interactions, 
respectively, and are active only when the two sites in- 
volved are occupied by different orbitals. The last term 
is AFM, with J3 = |t2[2/(c/2 + Jh) + l/(C/2 + 4Jff)], 
and is non-zero only when the two sites have the same 
orbital occupancy. Parameters that play a role in the 
Hamihonian (jTJ are: t = t„ ~ 0.32 eV, Jh ^ 0.64 eV 
and U2 ^ 4.1 eV.20'21 Thus tj = Jh/U2 ~ 0.15 < 1 and, 
just in order to present the results in a more transparent 
form, we expand the exchange energies around 77 = 0. 
We get Ji ~ J{l + r]), J2 ~ J{l-v) and J3 ~ 4:J{l-2r]) 
where J = /U2 — 25 meV represents the overall energy 
scale. 

The realization of the spin-orbital model ^ on the 
pyrochlore lattice has an immediate and important con- 
sequence, namely, that only some bonds can contribute 
to the energy, depending on their orbital configuration. 
Indeed, every bond ij in the a(3 plane has zero energy 
gain unless at least one of the two sites i and j has an 
occupied orbital of a/3 kind. The strength, as well as the 
sign, of spin-exchange energy associated with two NN 
sites i and j depends only on their orbital occupations 
and the direction of the ij bond. We can thus classify 
tetrahedral bonds in four types: 

(i) 60 J when both ions at sites i and j of the generic 
a/3-plane have a/3 orbital occupancy (Fig. la). It is char- 
acterized by a Hamiltonian with strong AFM exchange 
~ J: 

Hb„^-J{l-2ri){\~AS,-Sj). (3) 

(ii) 5i, if the two sites of bond ij in a/3-plane are oc- 
cupied by one 07 and one a/3 orbitals, 7 ^ /3(Fig. lb). 
This bond has a weak FM exchange ~ r]J and the corre- 
sponding Hamiltonian is: 

^ffci =-J(l + 7?/2 + 2r?5, -^O- (4) 

(iii) 62, with the two sites of bond ij in a/3-plane oc- 
cupied by one a7 and one (3^ orbitals. In this case there 
is no energy contribution: if^^ = 0. 

(iv) &3, if both sites of bond ij in a/3-plane are occu- 
pied by two a7 (/?7) orbitals. Also this bond does not 
contribute to the energy: iJbg = 0. 
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FIG. 1: Orbital arrangements on a bond: (a) bond 60 with 
strong AFM coupling ~ J, (b) bond 61 with weak FM cou- 
pling ~ 77J. 

Thus, only 69- and 5i-bonds have a binding energy and 
for any given orbital configuration on the pyrochlore lat- 
tice the Hamiltonian is just a linear combination of 
Hbo and Hb^. Moreover, the explicit form of Eqs. © 
and Q implies that, in the limit 77 ^ 0, a set of exact 
eigenstates of can be constructed: in fact, orbital in- 
teractions are already diagonal, and, in this limit, the 
only Heisenberg term that survives is that of Hba- This 
latter, by definition, is active only along the direction 
connecting the two sites with equally filled a/3 orbitals. 
This implies that the only "infinite" configuration of in- 
teracting spins is obtained when the spin-interactions act 
along the same direction, i.e., the AFM Heisenberg chain, 
which is exactly soluble i2S 

In the next subsection, starting from an isolated tetra- 
hedron we discuss possible orbital arrangements on the 
pyrochlore lattice and provide the solution of correspond- 
ing spin Hamiltonian. 



B. Possible phases on a pyrochlore lattice and 
their energetics 

It is possible to show that only three energetically in- 
equivalent tetrahedrons can be singled out of the four 
hn bonds introduced abovepiS, Topologically, they can be 
characterized by the number of strongest 60-bonds and 
classified in three following types [pictorially shown in 
Fig. 12]: A- type tetrahedra, with two 5o bonds, B-type 
tetrahedra with one 69 bond and C tetrahedra, with no 

60- bonds. Notice that, whatever the orbital configura- 
tion is, every single tetrahedron must obey the constraint 
2nf,Q -|- ribj^ = 4. This implies that A tetrahedra have no 

61- bonds, i3-type tetrahedra are characterized by 2 hi 
bonds and, for C tetrahedra, there must be 4 61-bonds. 
We can further classify S-type tetrahedra, according to 
the choice of the two orbitals on the tetrahedron bond 
opposite to the 60-bond, into: Bi, B2 and (see Fig. 
|2Jl where, respectively, the bond opposite to 5o one is of 
61, 62 and 63 type. We stress again that the number 
of 62 and 63 bonds is not relevant as far as the Hamil- 
tonian is considered: they acquire importance only 
when magnetoelastic correlations among bonds will be 
introduced, in the next section, as they allow to distin- 
guish topologically the various configurations in the py- 
rochlore lattice. A, B, and C tetrahedra are the bricks 
that allow to build the orbital pattern throughout the 
whole pyrochlore lattice. Because of the Ising-form of or- 
bital interactions, in the following we can focus simply on 
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FIG. 2: Orbital and bond arrangements on a tetrahedron for 
cases A, -61,2,3 and C discussed in the text. 



these three cases, relying on the fact that configurations 
with a linear superposition of orbitals on each site must 
have a higher energy. We shall do only one exception 
to study a case with a particular physical meaning, i.e., 
that of a "cubic" symmetry, where each site is occupied 
by a linear superposition with equal weights of the three 
orbitals -^[\xy) + \xz) + \yz)]. 

If we now try to cover the whole pyrochlore lattice 
with these tetrahedra, we can have the following global 
patterns: 

(I) Heisenberg chains: If all tetrahedra of pyrochlore 
lattice are characterized by two 60 bonds (i.e., all sites are 
occupied by the same orbital), then the effective Hamil- 
tonian ^ can be mapped into a set of one dimensional 
decoupled AFM Heisenberg chains that lie within the 
plane of the chosen orbital. The only interactions are 
due to feo-bonds and are described by the spin Hamilto- 
nian |(2Jl. Thus, the ground-state energy per site can be 
evaluated exactly by using the results for an Heisenberg 
chain^ that give Ea = -2.77(1 - 2tj)J. 

(II) Dimer phase B: This state is made of only B- 
type tetrahedra with one strong 60-bond, and two inter- 
mediate &i-bonds. As all three Bi tetrahedra are ener- 
getically equivalent all possible coverings of pyrochlore 
lattice by Bi tetrahedra have the same energy. When 
pyrochlore lattice is covered by -B^-types tetrahedra (two 
possible coverings are shown in Fig. |3J), then each spin 
is engaged in one strong AFM bo bond and two weak 
FM &i-bonds. Such coverings form a degenerate mani- 
fold and the corresponding energy can be calculated as 
follows. In the limit rj 0, the spin-only Hamiltonian 
can be solved exactly, as it can be decomposed into a 
sum of spin-uncoupled 60 bonds. In this case the en- 
ergy minimum is reached when the Heisenberg term of 
the 60-bond is the lowest, i.e., for a pure quantum spin- 
singlet {Si ■ Sj — —3/4). Remarkably, such spin-singlet 
(dimer) states, in the limit ry — s- 0, are also exact eigen- 
states of the full Hamiltonian (Q. As 77 <C 1, the dimer 
state is stable against the weak FM interdimer interac- 
tion. In this case the magnetic contribution along the 
FM &i-bond is zero {{Si ■ Sj) = for i and j belonging 



to different dimers) and we are led to an energy per site 
given by: Eg = Ebj2 + Eb,=-{3- ^77) J. Here Eb„^,^ 
is the energy of the bond 60(1). 

(HI) FM order: In this case all tetrahedra are of 
C-type with four interacting 61 bonds and two 63 (or 
one 62 and one 63) noninteracting bonds. Thus, all non 
zero spin-exchanges are ferromagnetic and given by Eq. 
Q . The ground state for this type of orbital ordering is 
thus ferromagnetic and has an energy per site given by 
Ec = 2Eb, =-2(1 + 77) J. 

(IV) Frustrated AFM: Even if the "Ising" form of or- 
bitals interaction implies that configurations with linear 
superposition of orbitals on each site must have a higher 
energy, we consider for completeness the case where each 
orbital is occupied by a linear superposition with equal 
weight of the three orbitals: -^{\xy) + \xz) + \yz)). The 
realization of this phase restores the full pyrochlore lat- 
tice symmetry and, after averaging Eq. over the 
orbital configurations on neighboring sites i and j, the 
system is described by the spin Heisenberg Hamilto- 
nian Hd/J = E»j[-5/9 + (4/9 - 1677/9)5, • Sj], whose 
ground-state energy per site is Ejj — —(1.89 — 0.897?) J. 
Here we borrowed the numerical evaluation for ground- 
state energy in a pyrochlore lattice from Ref. [23l| : 
{^/N)^--Si-Sj ~ —0.5. Such a state is highly frustrated 
and its ground state is a spin liquid»£ It seems worthwhile 
to note that the energy in this phase is higher than that 
of the other phases, as it does not exploit at all the po- 
tential energy-gain contained in the orbital orderingi^ 

(V) Mixed (dimer) phase AC: It is of course possi- 
ble to cover the pyrochlore lattice also by means of mixed 
configurations of tetrahedra. There are in principle infi- 
nite possibilities in this sense, but those minimizing the 
energy must contain, in average, at least one singlet bond 
every tetrahedron, due to the big energy gain related to 
the spin-singlet state in Hb,,. While it is not possible to 
have, in average, more than one spin-singlet bond every 
tetrahedron^ one can fill the whole pyrochlore lattice 
by means of a mixed configuration made of alternated 
tetrahedra with two and zero 60-bonds {A- and C- type 
tetrahedra, respectively). This configuration is shown in 
Fig2] For this type of orbital arrangement one can also 
construct an exact ground state of corresponding spin 
Hamiltonian, in the limit of vanishing Hund's coupling. 
The ground state is given by dimer phase where there are 
two spin-singlets, located opposite each other, on strong 
60 bonds of A-iype tetrahedra while C-type tetrahedra 
has no singlets. As in average, — 1 and — 2, this 
configuration is degenerate with Dimer phase B, as far 
as only dda overlap is considered. 



C. Ground State Manifold 

Energies of the possible phases on the pyrochlore lat- 
tice can be discussed in terms of r/, the only free parame- 
ter available. For ry = the lowest ground-state energy is 
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FIG. 3: The ground state coverings of the unit cubic cell 
through dimers. Locations of singlets are represented by thick 
links. Different numbers correspond to inequivalent tetrahe- 
dra. (a) The helical dimerization pattern (indicated by ar- 
rows) is formed by alternating short bo and long 63 bonds, 
Dimer phase B3. (b) One of the possible coverings of the 
cubic unit cell by B1B2 tetrahedra. 



that of spin-singlet degenerate manifold. With increasing 
rj we find only one phase transition at r^c = 2/11 ~ 0.18, 
from spin-singlct manifold to a FM phase. As 77c is above 
our estimated value of 77 ~ 0.15, we can conclude that 
the ground state manifold is given by spin-singlet phases 
spanned by degenerate dimer states of B or AC-types. 
This degenerate manifold is characterized by a static pat- 
tern of spin-singlets (dimers) throughout the whole py- 
rochlore lattice and, thus, is different from RVB state. 
Each dimer covering is frozen in an exact eigenstate of 
the Hamiltonian for 77 = 0. For finite 77 the different 
dimer patterns are not connected by the Hamiltonian: 
the bond corresponding to the dimer in each tetrahedron 
is fixed, being determined by orbital pattern and orbital 
degrees of freedom are static variables. Thus a tunneling 



FIG. 4: The coverings of the unit cubic cell through the tetra- 
hedra of A-type with two singlets and C-type with no singlets. 
Tetrahedra are labeled as in Fig. |S] 



between different dimer states can not take place. 

In this spin-singlet manifold the original spin degener- 
acy is removed. However there is still a remaining de- 
generacy to be lifted. The degeneracy of B-manifold is 
related to the freedom in the choice of the two orbitals on 
the tetrahedron bond opposite to the one of the singlet. 
Different choices of these orbitals give rise to inequivalent 
dimer covering patterns of the pyrochlore lattice with one 
dimer per tetrahedron (see Fig. This degeneracy is 
given by the number of such dimer coverings and the 
corresponding number of states can be estimated as fol- 
lows. In the S-manifold there is one singlet per tetrahe- 
dron and each spin can be engaged in only one singlet. 
Moreover when a singlet is located on a given tetrahe- 
dron then each neighboring tetrahedron is left with only 
three possible choices for a singlet location. Thus the 
number of such coverings grows with the system size as 

Ub ~ 3^^ = VS^. Here Nt = N/2 is the number of 
tetrahedra and we have ignored the contributions coming 
from closed loops (hexagons) on the pyrochlore lattice. 

The asymptotical determination of AC-manifold de- 
generacy is a more difficult task. However, based on the 
simple arguments we can estimate its lower and upper 
limit correctly. First, we divide the pyrochlore struc- 
ture in two sublattices formed by two differently ori- 
ented tetrahedra: one sublattice is composed by ^-type 
tetrahedron with two singlets located opposite each other 
and the other by C-type tetrahedron with no singlets. 
Then consider an A tetrahedron with two singlets in a(3 
plane: it can be directly verified that the four NN tetra- 
hedra of A sublattice, connected to it by straight lines 
in ck/3 plane, will have only two possible orbital choices. 
The remaining eight NN tetrahedra of the A sublattice 
have, instead, either two or three choices to locate the 
two spin-singlets. These constraints are dictated by the 
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condition that intermediate tetrahedron must be of C- 
kind. Thus, the degeneracy can be estimated as foUows: 
2^T < TV'ac < 3^^, where = N/A is the number of 
tetrahedra in one sublattice. We conclude that the de- 
generacy of this manifold is still macroscopic, but smaller 
than that of B one. 

Thus, even though the formation of spin-singlets re- 
moves the spin degeneracy of the pyrochlore lattice, there 
is still a macroscopic degeneracy to be lifted. The main 
question is whether this degeneracy can be removed by 
extending our effective Hamiltonian to the previously ne- 
glected ddiT and dd5 overlaps. When these processes are 
considered, 62 and 63 bonds acquire a different bonding 
energy, and, as number of these bonds is not the same 
for Bi and AC manifolds,-^ the degeneracy between these 
two phases can be lifted. However, this effect is smaller 
than the one induced by magnetoelastic coupling (see 

Sec. Ill), as Jdd-K = -jj^ — 2 meV, while magnetoelastic 
energy gain per ion is about 6.5 meV. It is then obvious 
to look for the degeneracy removal first in terms of this 
" spin- Teller" interaction, as done below. 

Moreover, the degeneracy within the S^-manifold can 
not be removed within the effective electronic model, not 
even introducing smaller NN hopping integrals. The rea- 
son is related to the fact that the energy gain depends 
only on the total number of bonds of each type (nho, 
nbi, ni,2, 7^63 ) in the unit cell, and, in order to fill the 
whole crystal with a periodicity not lower than the one 
of the primitive cubic cell, the average number of bonds 
per tetrahedron should be the same, whichever of 
the three tetrahedral configurations is taken {Bi, with 
n(,2=l and nb^=2; or B2, with nf,2=3 and nh3=0; or B3, 
with nf,2=2 and 71^3 = 1). This number is given by nb^ — l, 
nbi=2, n;,2=2, n;,3=l, and it corresponds to the value of 
B3 case, that is the only one that allows to cover the 
whole cubic cell without mixing to other configurations 
(see Fig. EJ. 

From the above discussion it follows that only correla- 
tions between bonds can lift this degeneracy. These cor- 
relations naturally appear if the magneto-elastic contri- 
bution to the energy is considered. The orbitally-driven 
modulations of the spin exchange interactions will distort 
the underlying lattice through the spin-Peierls mecha- 
nism and different distorted patterns will pay a different 
elastic energy. The three degenerate phases which will 
be discussed in the next section are those with the unit 
cell filled by all B3 tetrahedra ("Ba-phase"), that with a 
mixture of Bi and B2 tetrahedra ("B1B2 -phase"), and 
that where the unit cell is filled by alternated A and C 
tetrahedra (" AC-phase" ) . 



III. EVALUATION OF THE 
MAGNETOELASTIC ENERGY 

In a spin-Peierls system the magnetic energy gain due 
to the spin-singlet pairs outweights the increase in elastic 
energy due to the dimerization of the regular array. We 



have shown with qualitative argumentg^^ that in our case 
this mechanism can select the triplet-T normal mode of 
the tetrahedron group, leading to a distortion with one 
short and one long bonds located at opposite edges and 
four undistortcd bonds (see Fig. O. Indeed, a reduction 
of the bond length increases the magnetic energy gain 
and therefore favors the shortening of the bond with the 
strongest superexchange, i.e., bo, where the singlet is lo- 
cated. 

Even though this picture is correct in its basic fea- 
tures, yet a quantitative description of the global physical 
mechanism leading to the tetragonal distortion requires a 
more careful analysis, that takes into account all elastic 
normal modes of the single tetrahedron, as well as the 
correlations of these normal modes within the unit cell. 
The aim of the present section is just to analyze such a 
global mechanism. 

The dependence of the energy gain AE on the magni- 
tude of distortion can be evaluated as a sum of magnetic 
A£'™ and elastic Ai?*^' contributions at each bond: 

AE = J2^E^ + ^E:^- (5) 

V 

The sum is restricted to the NN sites. 

All bonds in the undistorted lattice have the same 
length and we represent their elastic energy assuming 
that all ions are connected one another through equal 
springs of constant k: 

AEtl = \k{5d,,f = l^^^f- (6) 

where the constant Co = kd'^ is the radial force 
constant^ and 5dij = c?.y — c?q is the deviation of the 
bond length from its value do in the undistorted cubic 
lattice. 

The magnetoelastic energy of the generic bond ij can 
be written as: 

AE^ = {J{da)+g5d,,%-S, (7) 

where g = ^^^^\d=do- The dependence of exchange con- 
stants on the distance is mediated by the hopping ma- 
trix element: exchange constants are proportional to i^, 
and t is inversely proportional to the fifth power of the 
distance. In Ref. |2y| it is possible to find a rough esti- 
mate of the proportionality constant, a: t = a/d^, with 

= jTldda—^- Here is a characteristic length, that 
for Ti is 1.08 A, while 'qdda = -16.2 and ^ = imeV- A^, 
giving the value a — —116.628 eV-A'^. It is immediately 
apparent that even a small reduction in the length of a 
given bond can lead to a relatively high magnetic-energy 
gain. 

In order to analyze the magnetoelastic behavior of the 
three degenerate phases of our spin-orbital model we first 
focus on a single tetrahedron, either of A, Bi or C type. 
For a more fluent reading, the detailed calculations are 



7 




FIG. 5: Triplet-T deformation mode from the irreducible 
representations of the tetrahedron group. This mode gener- 
ates a distortion of the tetrahedron, with short and long bond 
located opposite to each other and four intermediate (undis- 
torted) bonds. 



reported in Appendix B, and here we just comment on 
the main points. 

In general, the global magnetoelastic Hamiltonian in- 
volves more than one normal mode for each tetrahedron: 
for example, for A tetrahedra the singlet and one doublet 
modes contribute to the magnetoelastic energy, and for 
B2 tetrahedra all six normal modes are involved, as re- 
ported in Appendix B. Moreover, even though the triplet 
T-mode associated to the singlet bond leads to a bigger 
energy gain for B3 than for Bi tetrahedra, when all avail- 
able normal modes are considered, the total energy gain 
gets different contributions from all normal modes, but 

in such a way that it is the same for all three Bi tetra- 

2 

hedra, and equal to —4.5^. What makes the difference 
among these configurations in a specific case like that of 
MgTi204 is that when these tetrahedra are considered 
in an infinite lattice, the single-tetrahedron degeneracy 
is lifted by the correlations among tetrahedra. In fact, 
in correspondence with their common energy minimum, 
B3, Bi, and B2 tetrahedra are all distorted differently, 
as can be deduced from the location of their minima. 
Thus, when these distortions are combined together in 
order to fill in exactly the global cubic cell, as in Fig. 
El the appearance of new constraints, due to the rela- 
tive match of all tetrahedra within the cell, forces the 
solution to another minimum. For example, in order to 
keep the global volume unaltered, all bond distortions 
due to spin-singlets must appear with reversed sign from 
one tetrahedron to the other, because to each expansion 
there must be a corresponding contraction. Notice that 
also the AC-phase is still degenerate with Bi configura- 
tions, as far as just magnetoelastic energy gain associated 

with single tetrahedra is considered, as, in average, it is 

2 2 

again (— 8— 1)/2-^ = —4.5^ (see Appendix B). Remark- 
ably, when inter-tetrahedra correlations are considered, 
elastic singlet modes contribute to AC-phase differently 
of Bi phase. In fact, it is worthwhile for the system to ex- 
pand C tetrahedra and contract A tetrahedra of an equal 
amount, as this leads to a net energy gain, because the 
singlet mode energy gain of A tetrahedra exceeds that of 
C ones. 

It is clear from these considerations that in order to get 
the global magnetoelastic minimum for the unit cubic cell 




FIG. 6: Bidimensional projection of the unit cell. Circles 
represent Ti-ions, with the same labels as in Fig. m Di- 
ameters are proportional to the z coordinate of the ion: 
z — 1,0.75,0.5,0.25, in fractional units. All bonds are ori- 
ented along one of xy, xz, yz bisectors, as indicated. Helical 
chains (h.c.) are represented in the direction of increasing z. 
Doubly cut bonds are 60 bonds and singly cut are 61 bonds. 



we cannot rely on the unconstrained minimization of the 
single tetrahedra shown in Appendix B. One alternative 
possibility is to calculate the magnetoelastic energy for 
the normal modes of the global cubic cell for the three 
degenerate configurations {Bj,, B1B2, and AC-phase), in 
such a way as to automatically consider all constraints 
among tetrahedra. Yet, this procedure is not complete, 
as normal modes corresponding, e.g., to the buckling of 
tetrahedra do not contribute to the magnetoelastic en- 
ergy. Such a contribution appears only when the normal 
mode involves an increase/decrease of some bond length, 
and, thus, it can be expressed as a linear combination 
of the normal modes of each single tetrahedron, with the 
constraint that the position of each equivalent Ti-ion over 
all nearest neighbor unit cells be the same (k = modes) . 

In order to fulfill this constraint, we impose that the 
total length over all the straight bonds connecting one 
ion with the equivalent one in the next cubic cells is 
a conserved quantity. This means that, for example, 
dHi~H3 + dm-Ki + dK4~K3 + dK3~H4 (see Fig. is 
a conserved quantity, and, thus, the sum of the four dis- 
placements SdH4-H3+SdH3-K4 + SdK4-K3 + SdK3-Hi = 

0. The same happens for all other 3 chains in the xy- 
planc, for the 4 of the xz-plane, and for the 4 in the 
j/z-plane. Due to Eqs. ljB4p . these twelve conditions be- 
come constraints for the normal variables of the different 
tetrahedra. The search for the minima of the magne- 
toelastic Hamiltonians (|B5IB9ll with the previous con- 
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straints leads to the following expression of the energy 

per unit cell: -18^ for Bg-phase; -17.2^ for B1B2- 

2 

phase; —14.4^ for AC-phase. Thus, the ground state is 

2 

given by -B3, with an energy per tetrahedron of —2.25-^, 
corresponding to the sum of the energies of the triplet 
t2 and doublet ei modes. The singlet and the doublet 
62 modes do not contribute to the energy, when they are 

constrained on the global unit cell. As now each ion is 

2 

shared by two tetrahedra, the energy per site is —1.125-^. 

However, there is an experimental fact that is not 
properly explained by this solution, and for which an 
extra discussion is required. The solution of the con- 
strained minimization gives an amplitude of ^ for triplet 
t2 mode and | for doublet Ci mode, thus a ratio 2:1 
in favor of the triplet mode. But the experimentally 
reported tetrahedral distortionsiS are, referring to Eqs. 
(|53l . and IHII): 6ri3 = -0.155 A, Sr^ = -0.006 A, 
which imply, when only t2 and ei normal modes are ac- 
tive, ei = -0.012 A, and t2 = -0.149 A, with a triplet- 
to-doublet ratio of more than 10. This indicates that, 
in the previously outlined model, some features of the 
real system are not taken into account, like the elastic 
potential coming from the oxygens, that acts in such a 
way as to reduce some distortions. Thus, the possibility 
to get better quantitative results for the experimental 
bond-lengths within our scheme should passes through 
a way to mimic this aspect by means of a stronger con- 
straint on the available modes before the minimization 
procedure. If we consider Fig. El we can see that the 
ion iJ4 can be connected to the equivalent ones in the 
neighboring edge-sharing unit cells through one of the 
three "linear" paths: HA - HI - LA - LI - iJ4, in xz- 
plane, or Hi - H2 - MA - M2 - HA, in yz-plane, or 
HA — H3 — KA — K3 — HA, in xy-plane. Face-sharing 
unit cells are connected by the paths along the heli- 
cal chains, namely: HA — H2 — L3 — LI — HA and 
HA - HI - A/3 - M2 - HA, for chains along z-axis, 
HA-H2-K1-K3-HA and HA- M2- Ml- HZ- HA, 
for chains along a:-axis, HA — HI — K2 — K2> — HA and 
HA — H2> — L2 — L\ — HA, for chains along y-axis. In 
principle the global HA — HA distance along, e.g., the z 
axis can be kept fixed even if the sum of the four bonds 
changes, as the bonds are not along a straight line, and 
the elongation of one bond can be combined with an ap- 
propriate rotation, due to a tetrahedron buckling. Yet, 
this buckling implies an extra elastic energy loss due to 
the relative movement with Mg and O ions. We can 
forbid these extra losses, by imposing another series of 
constraints that express the global conservation of the 
total length also along the "helical" chains, in this way 
freezing some of the normal modes that were previously 
allowed. 

When the minimization is performed with this new 

constraint, the results for the energy per unit cell are: 

2 2 2 

— 16^ for i?3-phase; —13^ for _Bi_B2-phase; — |^ for 

the AC-phase. Again, Sa-phase is the lowest, but this 

time the amplitude of ei mode is zero and the minimum 



energy is reached for a pure triplet mode, of amplitude 
Notice that with this constraint the ^C-phase gains 
very few, and also the relative energy difference between 
-Si-S2-phase and B3 increases. 

Finally, let us stress again that the main result is 
that when correlations among magnetoelastic energy are 
taken into account, i?3 phase is always stabilized with re- 
spect to the two configurationally degenerate B1B2 and 
AC ones. Of the two criteria used for preserving the cell 
volume, the second seems more in keeping with experi- 
mental data. This suggests that oxygens play a role in 
the elastic distortions of the spinel cell, which is more 
rigid towards bucklings than the bare pyrochlore lattice. 
In any case it is found that the elongation of the non- 
interacting 63 bond and the corresponding stretching of 
60 bond throughout the whole cell is always energetically 
more favorable, as found experimentally. 

Given the previous d-dependence of the hopping am- 
plitude t, it is possible to estimate the gain energy per 
ion associated with this magnetoelastic distortion. First, 
g = -|^U=do = 10^. Using Eq. (El, ^^13 = -0.155 
A= with fc = ^. This implies that Co = ^§j^ ^ 

9.7 eV, if we take do — 3.008 A. Finally, the estimate of 
the average magnetoelastic energy per site can be easily 
found: there are 16 ions in one unit cell, and thus the 

energy per site is E^e = — Co ~ ^'^ meV. 

IV. EFFECT OF MAGNETIC FIELD 

As we have shown above the ground state of the sys- 
tem is non- magnetic and spins are paired in singlet states. 
The spin degrees of freedom are thus gapped and the gap 
of triplet excitations is given by the singlet binding en- 
ergy As — AJ[\ — 2rf\ (see Eq. Q). When an external 
magnetic filed is applied, the energy of triplet excitations 
decreases, due to the gain of Zeeman energy, while the 
singlet state does not experience the applied field. There- 
fore, one would expect no change in the symmetry of the 
ground state up to fields q^bHc = at which singlet 
triplet gap closes [g is the gyrotropic factor and the 
Bohr magneton). At H = He a second order transition 
can in principle take place from non-magnetic to a mag- 
netically ordered state, driven by condensation of lowest 
energy triplets at some ordering wave- vector Q at which 
triplet spectrum has a minimum. 

However, the presence of an orbital degeneracy can 
modify this convectional picture. As shown in Sec. II, 
the first excited state above Dimer phase is the FM one, 
characterized by C orbital pattern. The energy of FM 
phase, contrary to the Dimer phase, is sensitive to the ap- 
plied magnetic field and decreases of AEh = —gH^BSf. 
The critical field He at which the energies of the two 
phases are equal is: 

gfiBHe = J[2-Ur]], (8) 
and one can easily verify that He < He- Therefore at 
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H = He there will be a first order transition from Dimer 
phase to FM state. This transition will be accompanied 
by a simultaneous rearrangement of orbital ordering and 
the abrupt closure of the spin gap, as the FM phase is 
gapless. We can estimate the order of magnitude of He, 
for rj ~ 0.15 and J = 25 meV, as hbHc — 4.4 meV, i.e., 
a critical field of He — 76 Tesla. In this estimate we did 
not consider the extra stabilization energy in favor of the 
singlet phase coming from the magnetoelastic distortion. 
Of course, one can read Eq. JSJ the other way round, in 
terms of rj, and deduce that, if a critical field is found at 
a lower value than 76 T, this implies that JhIU2 is closer 
to the critical value of 0.18 than what estimated here. 



V. DETECTION OF ORBITAL ORDERING 
THROUGH RXS. 

As already outlined in the previous two sections, or- 
bital ordering in MgTi204 is mainly dictated by superex- 
change interactions. This implies that the orbital orien- 
tation at each site does not strictly follow the symme- 
try imposed by the local crystal field, as happens when 
the effect is purely determined by Jahn- Teller distortions. 
Therefore, we believe that there is the possibility to ex- 
ploit the local symmetry differences in the helical d^z- 
dyz orbital pattern by means of resonant x-ray scattering 
(RXS), where the local transition amplitudes are added 
with a phase factor that can compensate the vanishing 
effect due to the global tetragonal symmetry. Conceptu- 
ally, this procedure is the analogue of the one used in the 
case of manganiteaSi that led to a series of results initially 
interpreted as a direct evidence of orbital ordering, and 
soon later recognizedi^ii^ii^ as mainly determined by the 
oxygen distortion around each Mn-ion. Yet, in this case 
Bragg- forbidden reflections (e.g., for LaMnOa^^) relate 
two sites with different orbital occupancy and local dis- 
tortion. Thus, at K edge, the signal at a given energy 
turns out to depend on the difference between, e.g., Px 
and py density of states around that energy^SS projected 
on Mn-ions. Such local anisotropics in the electronic 4p 
density of states can be a consequence of oxygen distor- 
tions (Jahn- Teller effect) or can be induced by the or- 
dering of the underlying 3d orbitals through the id — Ap 
Coulomb repulsion Udv (orbital ordering effect) ii4 The re- 
sults of Rcfs. 13 14 15 independently show that the effect 
of the Jahn- Teller mechanism is much stronger, about a 
factor 10 in amplitude, than the effect of orbital ordering. 
This led to the conclusion that 00 cannot be "directly" 
probed by means of RXS at K edges, where the strength 
of the ligand field overwhelms the effect induced by U dp- 
Hue to this result, the search for experimental evidence 
of orbital ordering moved to L edge RXS'^*' where one is 
sensitive directly to 2>d orbitals. Yet, this kind of spec- 
troscopy has the serious drawback that, in order to have 
non-imaginary Bragg angles, it is necessarily limited to 
crystals with very big unit cells, and this is not feasible 
for MgTi204. In fact, Ti L2,3 edges are at about 453.8 



and 460.2 eV, which corresponds to a photon wavelength 
A ~ 27 A, that does not allow Bragg law sinOs — \/2d 
to be verified, not even for (001) reflection. 

In spite of all this, we believe that it is still possible to 
detect the proposed 00 at K edge. The key-point to go 
beyond the conclusions of Refs. p^3ll#15l | lies in the fact 
that their results are strongly related to symmetry and 
distortion of the system as well as the kind of reflection 
under study. In particular, while the signal induced by 
the Coulomb repulsion Udp is of the order of ~ 200 300 
meV, independently of the particular crystal structure;^ 
the influence of the oxygen distortion is much lower in 
the directions of t2g orbitals than in eg ones. It might 
indeed seem conceivable that a reduced Jahn- Teller dis- 
tortion can give rise to a sizable interference with the 
00-induced effect in the RXS intensity. Moreover, there 
is the possibility that the rotation of the exchange vector 
and/or polarizations in the RXS experiment can make 
the effect more or less pronounced, according to the ex- 
perimental conditions. We have thus performed a nu- 
merical simulation and found that it is indeed possible 
to experimentally reveal the presence of the 00 through 
the comparison of different RXS signals. In the following 
we illustrate the details of our calculations, in order to 
clarify the previous theoretical speculations. 

The transition process for RXS, governed by the Fermi 
Golden rule, depends on the state overlap through the 
polarized electric field of the incoming {i) or outgoing 
(o) photon. If we consider the multipole expansion up to 
the electric quadrupole contribution, we get: 



(9) 



Here if^g and ipn are ground and intermediate state 
wave functions, respectively, is the polarization of 
the incoming (outgoing) photon and fc'^"-* its wave vec- 
tor. Around an absorption edge, Af„g is highly energy 
and angular dependent and this provides the sensitivity 
to the electronic structure around the atom. In RXS 
the global process of photon absorption, virtual emis- 
sion of the photoelectron, and subsequent decay with re- 
emission of a photon, is coherent, thus giving rise to the 
usual Bragg diffraction condition. The outgoing photon 
can have different polarization and wave vector compared 
to the incoming one and the atomic anomalous scattering 
factor (ASF) reads^: 



hu) huo 

ng 



Here huj is the photon energy, nie the electron mass. Eg 
and En are the ground and intermediate state energies 
and r„ is the broadening of the transition. The sum over 
the intermediate states starts from the Fermi energy;^ 

In the tetragonal phase of MgTi2 04 the s pac e group 
is the chiral group P4i2i2, No. 92 of Ref. 1^ (or its 



10 



"mirror" -related P432i2, No. 96), with 8 ions per unit 
cell, and the structure factor is: 

^(Q)=^e''3-^~^7, (11) 
i=i 

where fj is the ASF of ion j. The atomic positions, 
their orbital filling and the symmetries relating the ions 
one another are schematized in Table I. Notice that the 
tetragonal cell is characterized by a 45-degree rotation in 
xy-plane, compared to the cubic cell of Fig. ^a.) . 



atom 


00 


position 


symmetry 


Th{m) 


dyz 


{u,v,w) 


E 


Ti2(L4) 


dyz 


{ — U, —V, ^ + w) 




Ti3(i/1) 




-v,^ + u,j + w) 




Ti4(L2) 




(i+w, i-u,|+w) 




Ti5(L3) 


dxz 




C2y 


TkiH4) 


dxz 


(i +u, i - w) 


C2x 


Th{H2) 


dyz 


{v,u, -w) 


C2XX 


Tig (LI) 


dyz 


{—V, —u, i — w) 


C2XX 



TABLE I. Orbital occupancy and tetragonal position of 
Ti-ions. Symmetry operations are referred to Tii . 

Here u = 0.9911, v = 0.2499, w = 0.8668 are the frac- 
tional coordinates of the atoms in unit of tetragonal axes. 
Notice that such a derivation allows also to deal with x- 
ray natural circular dichroism (XNCD), as the absorp- 
tion cross section corresponds to the imaginary part of 
the forward scattering amplitude (i.e., when (5 = 0). In- 
deed a signal in the dipole-quadrupole interference chan- 
nel could be expected, as the space symmetry group of 
MgTi204 is a chiral one, and there could be the possi- 
bility that it is affected by 00. Unfortunately, our nu- 
merical simulations with the finite difference method,^* 
in the same range of parameters as shown below for RXS, 
demonstrates that such a dichroic signal is too low (less 
than 0.05 % of the absorption) and that the changes in- 
duced by 00 are practically negligeable. 

For this reason in the following we shall analyze just 
K edge RXS, and in particular all reflections (00/), that 
are Bragg-forbidden unless I = An. Using the symmetries 
of Table I, the structure factor can be expressed in the 
following form: 

^(ooo = (i + (-)'C2.)(i + (^)'c'4+) 

Note that even if /i is a scalar, it is expressed as a 
scalar product of two tensors, and, in Eq. H12|l . ro- 
tation operators are understood as applied to one of 
them. We can limit to the dipole approximation, where 
we just need to consider rank-2 cartesian tensors, Tq/j, 



whose irreducible representation in SO (3) are a scalar, a 
pseudovector and a traceless symmetric tensor: the first 
never contributes to Bragg-forbidden reflections (which 
are scalar- forbidden), while the second is proportional 
to a magnetic moment and has no influence on non- 
magnetic systems like MgTi2 04. The third tensor has 
5 independent components, that can be labeled accord- 
ing to the usual second-rank real spherical harmonics as: 
Dxy, Dxz, Dyz, Dx2_y-2, and D^z^-r^- They are a mea- 
sure of the anisotropy of p-density of states projected on 
the resonant ion^ For example, D^2_y2 measures the 
difference in the density of states in x-direction, px, and 
in y-direction, py. In order to determine the implica- 
tions of Eq. (|12(l on the various reflections, we need 
to know how the symmetry operators act on our sym- 
metric tensors. The indices of each tensor change ac- 
cording to the following rules: C'2z{x, y, z) — (— —y, z); 
C2y{x,y,z) = {-x,y,-z) and Ctz{x,y,z) = {y,-x,z). 
Similar results are obtained reminding that the ASF in 
our approximation is a second-rank irreducible tensor 
and Ca./i'^ = (-)"V^?\ C+zfin^ = i^f!-^\ and 
C2yfrn = (-)'"/i™- As a conscqucuce of Eq. ^ 
we get that all Bragg-forbidden reflections of the kind 
(0,0,4n-|-l) and (0,0,4n-|-3) are sensitive to the complex 
mixture of density of states: Dxz ± iDyz oc Y^^. For this 
reason, due to the cylindrical symmetry, a constant az- 
imuthal scan is expected as well as a minor dependence 
on the orbital anisotropics of xz and yz kind. On the con- 
trary, Bragg-forbidden reflections of the kind (0,0,4n-t-2), 
if we neglect cos47rw ~ 0.1 <C sinAnw ~ 1, are just pro- 
portional to Dxy. thus, they must show a non-constant 
azimuthal scan, and, indirectly, manifest some properties 
related to the 00, through the depletion of the filled dxy 
orbitals. 

If we express the polarization dependence in the vari- 
ous channels, referring to the cubic frame of Fig. 1, we 
get the following results. Both (001) and (003) reflec- 
tions couple to ex^z and ej,e^, and thus they are different 
from zero only in the crTr-channel. As their amplitude 
is proportional to (sin0 -|- zcos0), they have a constant 
azimuthal scan. Instead, (002) reflection is detectable 
in all aa, air and tttt channels, with an intensity that 
scales, respectively, as sin^(20)_D^y, sin^ 9b cos"^ {2(j)) D'^y , 
sin'^ 9b sm^{2(j))Dly, where 9b is the Bragg angle. As 
sin^ 9b — 0.1, signals in TTTr-channcl are 1/10 smaller than 
those in air and 1/100 smaller than those in (Tcr-channel. 
Also their azimuthal scans are out of phase: that of (ttt 
channel has its maximum value when <^ = 0, i.e., in the 
direction of the nearest neighbor oxygens, while those of 
aa and tttt channels take their maximum value along the 
Ti-Ti chains in the xy-plane, where an signal is zero. 

In Fig. 13 we show the results of our numerical simula- 
tions for (OOZ) reflections, performed with the finite dif- 
ference method option of the FDMNES programiiS As a 
input file we used the refined positions of Ref. for 
Mg, Ti and O ions. Because of the increasing CPU-time 
consuming, we performed our simulation, with and with- 
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FIG. 7; RXS intensities for (00/) reflections = 1, 2, 3) with- 
out 00 (a) and witii OO (b). Angles refer to incoming g po- 
larization: (j) = Q° coincides with the cubic x-axis and (j) = 45° 
with the tetragonal x-axis. The cluster radius is 3.3 A. The 
reference energy is the Fermi level. 



out orbital ordering, just for a cluster of radius 2.5 A, 
and one of radius 3.3 A, the first including the oxygen 
octahedron and the second also all 6 nearest neighbors 
Ti-ions around the resonant atom. All assertions stated 
before, regarding the relative intensities and polarization 
dependence, have been confirmed numerically. 

The main result of our numerical simulations lies in the 
comparison of Fig. Ufa) and Fig. |3[b), that shows how to 
reveal the presence of 00 through a simple experiment. 
In fact, the 00 case is characterized by a big increase 
in intensity of (002) refiection, of a factor of 1.7, while 
(003) is left unaltered and (001) increases of no more than 
1.3. Thus, an experimental setup that measures these 
signals keeping track of the relative intensity is sufficient 
to determine whether the orbital order is present or not: 
if the three reflections have almost the same intensity, 
there is no 00, while if the (002) is about twice than 
the others, this is the signature of the 00 of the kind 
predicted by our theory. 

In order to understand the reasons behind this en- 
hancement in intensity, it is useful to study K edge ab- 
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FIG. 8: Absorption spectra for a TiOe cluster in MgTi2 04. 
Directions of incident polarizations are expressed in the 
tetragonal frame. Energy zero is at the Fermi level. For each 
direction the comparison is made between the OO case and 
that without OO. 



sorption, with and without 00, in the fictitious system 
where the only absorption center is Tii-ion. In fact, in 
this way we probe directly the p density of states pro- 
jected on Til-ion, i.e., /i, the square of which is the RXS 
intensity shown in Fig. [7| Thus, linear dichroism, e.g., 
along (100) and (010) tetragonal directions is sensitive 
to Dxy In Fig. |Slwe report the absorption for a single 
Ti-ion surrounded by the oxygen octahedron. Oxygens 
in xy-plane almost lie along cubic x and y axes, slightly 
squeezed towards tetragonal x' axis. Labels in Fig. |S1 
refer to the tetragonal system, so that (100) direction 
corresponds to the incoming polarization along the line 
connecting the nearest neighbor Ti-ions in xy-plane (45° 
from the x-axis of Fig. It is evident that in all figures, 
except Fig. |HIa), the presence of 00 does not change the 
signal and the dichroism is entirely due to the ligand dis- 
tortion. In particular. Fig. |SJc) shows the dichroism 
in the yz cubic plane, while Fig. Eld) is taken in xz 
plane and Fig. Elb) is a measure of Px — Py, towards the 
oxygens. On the contrary. Fig. |SIa) shows, in correspon- 
dence to its flex, an offset between the two curves varying 
from about 0.6 eV, in the absence of 00, to 0.9 eV, when 
00 is present. 

Finally, we want to comment about the influence of the 
cluster radius on the results. A series of numerical sim- 
ulations of the local absorption (like those of Fig. |SJ) has 
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been performed with radii 2.5, 3.3 and 4.0 A, by means of 
the muffin tin program. The results show that the maxi- 
mum dichroism between the (100) and (010) absorptions 
due to Jahn- Teller corresponds to the 3.3 radius, where, 
compared to the 2.5 cluster, two Ti-ions are added along 
the (100) tetragonal direction, whereas no Ti ions are 
present along the (010). The simulations with the 4.0 A 
cluster shows again a reduction of dichroism compared 
to 3.3 A, due to the cell composition. Thus, the results 
shown in Fig. [7|for the cluster radius of 3.3 A correspond 
to the theoretical lowest ratio between the 00 signal and 
JT one, and, thus, the effect of 00 is expected to be not 
lower than that shown in Fig. |3[b). 
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With these premises, the effective spin-orbital Hamil- 
tonian can be written as: 

Hcs = -jj^^Y. ■ + 3/4] Of^ 

+ Tf;h^T.[S.-S,-l/i]Of^ (A2) 

+ J2[S.- S, - 1/4] {D^O^/ + D^O^/) 

(ij) 

where the sum is restricted to the NN sites of the py- 
rochlore lattice. The first term in HeS describes ferro- 
magnetic interactions, while second and the third terms 
arc antifcrromagnetic. The explicit form of orbital con- 
tributions depends on the kind of intermediate excited 
states involved in the supcrexchangc process: when the 
intermediate states have just singly-occupied orbitals 
{of A, their expression is: 



APPENDIX A: THE "FULL" EFFECTIVE 
HAMILTONIAN. 

In a previous Letteni^ we have treated the effective 
Hamiltonian a la Kugel-KhomskiiAi with dda contribu- 
tions only. In this appendix, for completeness, we take 
into account also the smaller terms proportional to the 
orbital overlaps of ddn and dd5 type. In this case the 
hopping matrix is no more diagonal as in Eq. [3 If we 
label planes and orbitals directions as j/z = 1, = 2 
and = 3, the matrix elements can be expressed in the 
form, given by Slater and Kostes^^: 

= Z/Adda + 1/Add5 
tf^^^^ = l/2dd^ + l/2ddS (Al) 

,ab _ n 

i^^^^^j = l/2ddTT~l/2ddS 

where a, b and c can be 1, 2, or 3. Lower indices denote 
the plane (a) where the bond ij lays, and upper ones label 
the orbitals. For example, t^jf^^-^ is the hopping matrix 
element from yz to xz orbitals along the bond ij in xy- 
plane. 

Then, we introduce on-site transition orbital operators 
Tj"'', which describe initial and final orbital states in the 
hopping process at site i. If there is no change in the or- 
bital state at site i, then operators are simply orbital 
projectors and can be written as: T/^ = — ^rf (1 — rf ), 
TP = irf (1 + rf ) and Tf = (1 + rf )(1 - rf ). Non- 
diagonal matrix elements of the transition orbital op- 
erator describe changes in the orbital state during the 
superexchange process: T/^ = i'^'i^'^'j^i = \t~ t~ , 



3 3 3 3 

a=l 6=1 c,d=l c=l 
3 3 

+ E E ^f^TfTf (A3) 

a,b=l c,d = l 

This orbital term is common to both FM and AFM 
bonds, whose difference resides in their energy denomi- 
nators, l/(C/2 - Jh) for the FM bond and 1/{U2 + Jh) 
for the AFM. 

The third term of the Hamiltonian describes the anti- 
ferromagnetic superexchange evolving through an inter- 
mediate excited state with double occupancy. The partic- 
ular form with two energy denominators {Di and D2) and 
two corresponding orbital terms is originated by the fact 
that orbitally doubly-occupied states are not eigenstates 
of the on-site multiband Hubbard Hamiltonian;^ due 
to interorbital double hopping. Exact eigenstates with 
single-orbital double occupancy are just linear combina- 
tions (with plus or minus sign) of two orbitally doubly- 
occupied ionic states, and have two different energies, 
U2 + Jh and U2 + 4 Jjj. If, in the superexchange process, 
the electron forms a double occupancy at some orbital 
on site j and comes back from the same orbital, then 
the energy denominator is Di = 3^^^^^^^ + jijjp^y If, 
on the contrary, the orbital state on site j changes dur- 
ing the hopping process, then the energy denominator is 
^2 = ~"S{uT+J) + 3(^2+47) ■ The corresponding orbital 
contributions to -ffoff are: 

3 

0^'= E ^^rTf (A4) 

a,6,c— 1 
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(A5) 



a,fc=l c.d=l 
c^b 



In the absence of ddn contributions, the hopping ma- 
trix connects only orbitals of the same type, a, which 
have nonzero overlap only in the corresponding a plane. 
Then the effective spin-orbital Hamiltonian <jA2|l simpli- 
fies to the one we have obtained in Ref. and which we 
deal with in Section II. 



APPENDIX B: MAGNETOELASTIC NORMAL 
MODES FOR A SINGLE TETRAHEDRON. 

The normal modes of a single tetrahedron are a well- 
known argument: yet, as we need their explicit expres- 
sions both here and in Sec. Ill, we prefer to review the 
basic results. Considering the small ionic displacements 
around the equilibrium positions, the coordinate of the 
4 ions, when the origin is at the center of the tetrahe- 
dron, read: ii = {do + xi,do + j/i, rfo + ^i), 12 = (—do + 
X2,-do + y2,do + Z2), ia = (-do + 0:3, do + 2/3, -c?o + 23), 
14 = (do -|- X4, —do + j/4, -dp + Z4). Assuming that the 
four ions are connected each other with six equal springs 
of constant k, there will be an associated potential en- 
ergy expressed as a bilinear form of the 12 displacements 
Xi,...,Z4. The diagonalization of this form leads to the 
12 normal modes of the tetrahedron, six of which are the 
zero-energy modes associated to the three global transla- 
tions and the three global rotation, that we may neglect 
in the following. 

The remaining six eigenvectors can be classified into 
a triplet T-mode, with eigenvalue kt = 4fc, a doublet E- 
mode, with eigenvalue ke = 2k and a singlet S-mode, 
with eigenvalue fcs = 8fc. They all correspond to a par- 
ticular deformation of the tetrahedron: those of T-modes 
are shown in Fig. \S\ those of the E-mode correspond to a 
local tetragonal or orthorhombic distortion (see, e.g., Ref. 
0), and that of the S-mode is the "breathing" mode, i.e., 
a volume expansion or contraction that does not change 
the tetrahedron shape. If we introduce the stretches of 
each bond, Sdij, in terms of the ionic displacement coor- 
dinates, it is then possible to express normal coordinates 
in terms of these latters. 



Sdi2 = 

Sd34 
ddi3 : 

Sd2i : 

(5di4 



1 

71 
1 

71 
1 

71 
1 

71 
1 

71 



(xi -X2 + yi- V2) 

(X4 - X3 + ^3 - y^) 
(Xi -XS.+ Zi- Z3) 
(X4 - X2+ Z2- Za) 



(yi - y4 + - z^) 



(Bl) 



and 



'5d23 ^ --j={y3 - y2 + Z2 - 2:3) 



-^{5di2 -t- (5r34 + (5di3 + 5d2i + Sdu + 5d2z) 
2V6 

i(5di3 + ^d24 - 5dii - 5d23) (B2) 



1 



2%/3 
1 



{25di2 + 25dzi — Sdi3 — 6d24 ~ Sdi4 — Sd23) 



yt2 



7*3 



((5di2 - ^d34) 



((5di3 - Sd: 



2i) 



((5di4 - (5d23) 



The inverted relations are: 



134 



Jl3 































).i 








1 






2' 


hl + 






1 








2' 


hl- 






1 








2' 


hl + 






1 








2' 


hi- 



(B3) 



7*3 



The spin-lattice interaction Hamiltonian becomes, 
then: 



\/6 



E: 



34 



Ei3 + _B24 + £'14 + £23) 



-gOei{Ei3 + E24 — Ei4 — E23) 
1 

gO e2{2Ei2 + 2i?34 — Ei3 — E24 — Ei4 — E23) 

(B4) 



2%/3 

gOti{Ei2 — E34) — g6t2{Ei3 — E24 
gOt3{Ei4 — E23) 



where Eij denotes the energy of the ij bond in J units 

and g = ^^^|d=doi with do the equilibrium distance. 

In this way it is straightforward to have the energy gain 
related to each mode of our spin-orbital model. Treating 
separately all the cases of topologically distinct tetrahe- 
dra, we get: 

1) A-tetrahedron with two spin singlets located at 
bonds di2 and d34. Superexchange energies are: E12 = 
E34 = -4 and Ei3 = E24 = £14 = £23 = 0. The global 
magnetoelastic Hamiltonian HmeA — Hsia + H^ia is: 



14 



HfneA — 



16 

71^ 



-keel. 



(B5) 



As each normal mode is independent of the other, it 
is possible to minimize with respect to the 2 variables 
separately and get the energy minima: Esm — ~^^~-< 
and i?e2m — The normal modes corresponding 



8i? 



to these minima are: x^sm — yj^: ^e-zm - 
As ks = 2kt = 'ike = 8fc, the global minimum is: Em — 

2) _Bi -tetrahedron with superexchange energies: £'14 — 
E24 = -1, Ei3 = -4 and E23 = E12 = E34 = 0. The 
global magnetoelastic Hamiltonian is: 
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1 

— / 
2 

1 



H„,eB2 = ^g0s + i;ks0l + 2geel + -keOli 



V3gee2 + ^kedl 



(B6) 



390t2 + ^kte^2 



9*3 + T;kt0l 



- V3gee2 + IkeOl, 



%2 + \ktel, (B8) 

The energy minima are: Esm. = -12|-, iJei™ = 

o 2 2 

i'fi2m = ~5f~j Eind i?t2m = ~8f7. The normal modes 

I2g 



corresponding to these minima are: 9s 



-'elm 



V3|-, and 6*; 



t2m 



-4f . The 



global minimum is still: Em — ~4.5^. 

Due to the energy gain in the triplet i2-iiiode, there is a 
shortening of the bond where the singlet is located, as it 
can be seen by substituting back the previous expressions 
for normal minima in Eqs. (|B3|l - IIB4p . 

It is worthwhile to note that even though all three B 
tetrahedra have the same magnetoelastic energy, the cor- 
responding bond distortions are different, and thus the 
three configurations match one another differently when 
forced in the global fee cell. 

5) Ci -tetrahedron (the one shown in Fig. |2Ie)) with 
superexchange energies: E12 — E34 — and E13 — 
E24 — Ei4 — E23 = —1. The global magnetoelastic 
Hamiltonian is: 



The energy minima are: £'s 



Ee2m = -ff;, Et2m = " ff; , and Et3m = -^r- The 
normal modes corresponding to these minima are: 9gm = 
2^ a. - „ ^ Tax, - -3f , and 



— —0-3- l-i 

— /c ' e2m 



The global minimum is: Em = —4.5— 



3) i?2-tetrahedron with superexchange energies: £^14 = 
E12 = -1, Ei3 = -4 and £;23 = £^24 = E34 = 0. The 
global magnetoelastic Hamiltonian is: 



HmeBS — 

V3 



\ksel + |g0ei + \keel. 



2 ,^e2 + \keel2 + geti + \hel^ 
+ ^get2 + \ktej2 + 9et3 + \ktOh (B7) 

2 q 2 

The energy minima are: Esm = -12f7, -Beim = -gf;, 



£. 



3r 

8fee 



Et, 



and £;tim = £t3m ^ 



The normal modes corresponding to these minima are: 

12g r\ 3g n V3g n 



and ^tim = Qtzm = ~k^- The global minimum is again: 

Em = -4.5^. 

4) i?3-tetrahedron with superexchange energies: E14 = 
£23 = -1, Ei3 = -4 and £24 = £^12 = £^34 = 0. The 
global magnetoelastic Hamiltonian is: 



sBl — 



12 



-k fP' 

tefel 



meCl 



2 



1eCe2 



(B9) 



The energy minima are: Esm = and Ee2m = 

2„2 

— The normal modes corresponding to these min- 
ima are: 0sm = ^^"^ ^^2™. = The global 
2 

minimum is: Em — 

6) C2 -tetrahedron (not shown in Fig. 01 • They are 
characterized by 4 61-, one 62- and one 63-bonds, with su- 
perexchange energies: E23 = E24 = and £13 = £14 = 
E12 = E34 = —1. The global magnetoelastic Hamilto- 
nian is: 



HmeC2 — 



:ks0i 



1 



: fee 6*62 



yt2 



7f3 



:kt0t3 



(BIO) 



The energy minima are: Esm = £e2m 



and Et2v 



E, 



9 

2 

-T^. The normal modes cor- 
responding to these minima are: 6'sm — 

0e2m = -Ati and 9t2m = 0fim ^ ~ — 



and 

yg^, J ci±iu i^£2m — "tiJn — fcj"' 

case of i?i tetrahedra, the partial contributions of each 

mode compensate and the global minimum is, as for Ci 

2 

case, Em^ 
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